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Abstract 

The wave-packet treatment of neutrino oscillation developed previously is extended 
to the case in which momentum distribution functions are taken to be a Gaussian form 
with both central values and dispersions depending on the mass eigenstates of the 
neutrinos. It is shown among other things that the velocity of the neutrino wave packets 
does not in general agree with what one would expect classically and that relativistic 
neutrinos emitted from pions nevertheless do follow, to a good approximation, the 
classical trajectory. 
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§1. Introduction 

In a previous publication, 0* we have developed a wave-packet treatment of neutrino oscil- 
lation & for each of the "equal-energy prescription" and "equal-momentum prescription" (see 
Sect. 3 below), and, by invoking relativistic kinematics as well, have derived the necessary 
conditions for oscillation to occur, which appear to have a form more well-defined and quan- 
titative than what have been noted before. El The wave packets corresponding to each of 
the mass eigenstates of neutrinos in the case of the equal-energy (-momentum) prescription 
have been constructed from momentum (energy) distribution functions of a step-function 
type having a common central value and a common dispersion. 

In the present note, we like to extend our wave-packet treatment of neutrino oscillation 
by introducing momentum distribution functions with both central values and dispersions 
depending on the mass eigenstates of the neutrinos and by replacing step-functions by Gaus- 
sian functions. Such extensions enable our treatment to have closer contact with actual 
experimental or observational situations, though our emphasis is still placed more or less on 
conceptual aspects of neutrino oscillation. We shall see among other things that the velocity 
of the neutrino wave packets does not in general agree with what one would expect classically 
and that relativistic neutrinos nevertheless do follow, to a good approximation, the classical 
trajectory. @* We shall also argue that the equal-energy (-momentum) prescription seems 
appropriate to approximately describe oscillation experiments which involve something cor- 
responding to time- (space-) integration or averaging. 



§2. Wave-packet treatment 

Let \u a ) (a = e,/j,,r, etc.) be neutrinos associated with electron, muon, r lepton, etc., 
which are mutually-orthogonal superpositions of the mass eigenstates \ui) having mass mi (i 



1, 2, 3, etc.): 



W, 



= (2-i; 



Suppose a neutrino of flavor a is born at t — and propagates towards the x-direction. 
Then, its state vector at x and t may be written asB-* & 

\v a (x,t)) = ^2U ai gi(x,t)\vi) , (2-2) 



9iM = \l^- r dpJMe^*-^ , (2-3) 



2tt 



2 



where fi(pi) is the amplitude for creation of \ui) with momentum p$ and energy Ei = Ei(pi) 

(2-4) 



pf + mf , normalized aa^l 

/oo 
dpi\fi(pi)\ 2 = l. 
-oo 

The probability to find a neutrino with flavor a' at x and i is calculated as 



P a _> a /(x,t) = \(v a >\u a (x,t))\ 2 



where 



Gij(x,t) = g i (x,t)g j (x,t)* 

1 p QQ p OO 

Z7T J —oo J —oo 

As fiijp), we take, for simplicity and definiteness, 



fiiPi 



1 r fe-rf) 2 ! 

■exp 



Expanding Ei(pi) around p^fcf) 



where 



.0/irO 



# = (dEt/dp^o = p°/E 
and performing integration over p i5 one obtains 



g i (x,t)=il> i {x-ftt)e«i<» l! > 



where 



1 . x l 

exp[ 



6°(x,t)=p°x-E°t } 



*' See Appendix A for some related remarks. 
**) If one expands Ei(j>i) as 



(2-5) 



(2-6) 
(2-7) 



(2- 



(2-9) 

(2-10) 
(2-11) 



Ek ~ £° + (dEi/dpi^ofa-p?) + {d 2 E l /dpi) p ^ p o{p l -p'i) 2 /2 , 

one may discuss spreading of the wave packets with time.uP This is however beyond the scope of the present 
investigation. 
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Ef = ^{pfj 2 + mf and Ax,, = 1/Ap { . Substituting Eqs.(2.9) ~ (2.11) into Eq.(2.6) and in 
turn into Eq.(2.5), one gets 

G l3 (x,t) = ^(x,t)e 10 ^ , (2-12) 



t) =i/)i(x- $t)^{x - $t) 



exp[ 



% Ax i Ax j 1 2{Axi) 2 2(Axj) 2 1 



exp[- 



\-KAxiAxj 2(Ax i ) 2 (Ax j ) 2 2((Ar,) 2 + (Ax 3 ) 2 y ' 

(2-13) 

eJ(M) = 0?(s,t)-^(a;,t) 



(p°-p>-(£°-£°)*, (2-14) 



an> 



P a -,a'(x, t)=J2 UlPl^uix, t) + 2 J] UaiUanUajU^j&ijix, t) COs6>°(x, t) , 

(2-15) 

where 

?0 _ /3?(4w) a + /? (4Pi) a 



(Z^) 2 + (Z^) 2 



(2-16) 



l/CAr*) 2 + l/(Ar,) 2 ' 

We shall hereafter refer to \Pij(x,t) and 6>^(x, t) as the wave-packet factor and the phase 
factor, respectively. 

Trajectory of the neutrino wave packets 

\Pij(x, t) as a function of x has a peak at 

x = /??•*, (2-17) 

and one has 



/ T (flp _ fl0\2j2 

3KHtJ ' ; \nAxiAxj Fl 2{(Axi) 2 + {Axjfy ' v ; 



= (t^ai) is taken to be a real orthogonal matrix here. 
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+ ho I EW77^ I TaZ \2\ 1* 



l E? + E° (E? + E$)((Axi) 2 + (Axj) 2 ) 

(2-19) 

Eq.(2.18) tells us that ^(/^i, t) would reduce by a factor of e or more compared to its value 
at t = unless 

- $$)\t < p{{Ax t y + {Axtf) . (2-20) 

This gives the coherence condition for neutrino oscillation. Similarly, &ij(x, t) regarded as 
a function of t has a peak at 



t = xjffi , (2-21) 



with f3fj given by 



A 0/ 



i.i 



mf{A Vl f + ^f{A Vj f 

®{A Pl y + (3»{Ap 3 y 

ff/iAxtf + ^/iAxrf 



(2-22) 



and one finds 



Eq.(2.23) gives another form of coherence conditions for neutrino oscillation: 



]x . 



(2-24) 



Space- and time-integration prescriptions 
Let us now calculate 

/oo 
dxG i:i {x,t) . (2-26) 
-oo 

Substituting Eqs.(2.12) ~ (2.14), and performing integration over x, one obtains 
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r (a° - m 2 

xexp[ 2{{Ax l y + {A Xj yy 

^-V^Ax^Ax,)^ 

x exp[ ~ 2((A Xl y + (Ax,y) ] ■ (2 ' 27) 

Note that the same phase factor as Eq.(2.19) appears here. In the right-hand side of 
Eq.(2.27), the second line gives the coherence condition (2.20), while the third line gives, 
when p® ^ p®, another condition for neutrino oscillation to occur. These conditions, in the 
case of Api (or Axj) being independent of i, read 

\($-$)\t/2<Ax = l/Ap<2/\p° i -p° j \ . (2-28) 

Calculating 



G i:j (x) = I dtGij(x,t) (2-29) 

in a similar way, one obtains 



Gij(x) 



2(Ax l /Pf)(Ax J //3^) pi&Ux>x/ ^ 
X\ ^((Ax i //3ff + (A Xj /flr 
(1/09 - i/fi) V 



x exp[ 



2((A Xi /p?y + (Axj/fly 



x exp[ w^im + (A^im 1 ' ( } 

where O^Ax, x//3^) is given by Eq.(2.24). A result similar to Eq.(2.30) was obtained before 
by Giunti, Kim and Lee& for the case of Api (or Ax{) being independent of i, who noted 
that, in the right-hand side of this equation, the second line gives the coherence condition 
Eq.(2.25), and the third line acts as another factor to suppress the oscillation probability. 



§3. Comparison with the conventional approaches 

In the conventional approach, it is supposed that \v a (x,t)) is given by 

K(x,t)) = Y,u*My 9 ^ , (3-i) 

i 

and, accordingly, P a ^ a i (x, t) is given by 

P a ^ a> (x, t) = J2 UMi + 2 E U ai U aH U aj U alj cos 6>° (x, t) , (3-2) 

i i<j 
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where 8®(x,t) and &%{x,t) are given by Eqs.(2.11) and (2.14). If one further assumes either 
Ef = E® (equal-energy prescription) or p° = p® (equal- momentum prescription), one has 

0%{x, t) = (p* - p])x = -K 2 - m 2 )x/(rf + p]) , (3-3) 

or 

= -(£° - Efj)t = -{m\ - m))t/{E« + £°) , (3-4) 

implying that P Q ,^ Q ,/(x, t) will oscillate as a function of x or t with wave length £ij = 2tt|(p° + 
p°)/(m 2 — m 2 )| or period = 27r|(i? l + E®)/{mj — m 2 )|. Neutrino oscillation is usually 
discussed on the basis of Eqs.(3.1) ~ (3.4) and it is often claimed that these two prescriptions 
give practically the same results for relativistic neutrinos.0 

With the plane-wave expression (3.1), nothing can be said about the trajectory of the 
neutrinos. If, however, one arbitrarily assumes that the neutrinos are on the trajectory given 
by 

x = fat , fa = (p° + P ;)/(£° + E*) , (3-5) 
one would then find expressions exactly same with Eqs.(3.3) and (3.4): 



»;,(■>■■■>■/ hj) = "(m 2 - m 2 )x/(rf +p°) 
OU^t, t) = -{m\ - m))t/{E« + E$) . (3-6) 

We shall refer to the trajectory described by Eq.(3.5) as "classical trajectory" & and to such 



an approach as the " center-of-mass velocity" prescription. 
A couple of comments are in order. 

(1) In calculating Gij(x), Eq.(2.29), if one substitutes Eq.(2.7) and performs integration 
over t first, one would have a 5-function with its argument given by E^ — Ej, which implies 



"equal energy" explicitly.Ei3 Therefore, one may in this sense regard the conventional equal- 



energy prescription reviewed here as a simplified version of the time-integration prescription 
considered in the preceding section and regard the second term in Eq.(2.24), in the case of 
Ef ^ Ej, as a correction term to Eq.(3.3). Similar comments apply to the equal-momentum 



*' In our previous paper, on arguing that these two prescriptions are better to be conceptually distin- 
guished from each other, we have, in order to distinguish from the oscillation length £jj, referred to 
as oscillation period. As for distinction between these two prescriptions, see also the arguments given by 
Lipkin.& 

**) In Ref.8, the authors considered a frame defined by the velocity (3^, i.e., the center-of-mass frame of 
the two components of the particle in question^ 

***) In our previous wave-packet treatment, IP we have imposed the condition of equal energy by hand, 
as usually done in the case of the plane-wave approach, and thereby derived the conditions for neutrino 
oscillation to occur, i.e., the conditions corresponding to the inequalities L2J28). These conditions have been 
well known and derived mostly from more or less intuitive arguments. W u 
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prescription mentioned here and the space-integration prescription given in the preceding 
section. 

(2) The trajectory of the neutrino wave packets given by Eqs.(2.17) and (2.16), or given 
by Eqs.(2.21) and (2.22), differs in general from the trajectory given by Eq.(3.5). {flfj) 
given by Eq.(2.16) (Eq.(2.22)) would coincide with /3fj given by Eq.(3.5), and the phase factor 
given by Eq.(2.19) (Eq.(2.24)) would coincide with that given by Eq.(3.4) (Eq.(3.3)), if the 
central values and the dispersions of the momentum distribution functions of the neutrinos 
satisfy 

A Pi /A Pj = ^Ef/E] ( fiApi/fiApj = JpWj ) • (3-7) 
We shall shortly see that these relations do not necessarily hold. 

§4. Neutrinos from pion decays 

Now let us consider the case in which the neutrinos in question are created as a result of 
the decay 7r + — > fi + + v a (a = //) (with a branching fraction of 100%), and suppose that a 7r + 
with momentum p^ and energy E n = E 7r (p n ) = \jp\ + m%, moving towards the ^-direction, 
emits at t = a v a towards the x-direction (and a with mass m M , towards the opposite 
direction) in its rest frame. The momentum pi and energy Ei of the Vi component of the v a 
are given by 

= lM + P*E*) = (E w p* + p„E*)/m n , 
Ei = ln (E* + (3 nP *) = (E W E* + p„p*)/m„ , (4-1) 

where /3 n = p^/E^ and 7^ = l/yjl — f3% = E n /m n are the velocity and 7-factor of the 7r + , 
and p* and E* are the momentum and energy of the V; L in the rest frame of the 7r + , 

p* = M*) 2 - ™ 2 

= \J[ml - (m^ + mi) 2 }[ml - (m M - m;) 2 ] /2m w . (4-2) 

Pi and Ei are functions of p n and one readily verifies 

dpi/dp n = Ei/E w , dEi/dp n = Pi/E n . (4-3) 

It is natural to suppose that the momentum distribution fi(pi) of z/, is determined by the 
momentum distribution f n {p n ) of the parent n + and to formulate this relation quantitatively 

as 

\fi(pi)\ 2 dp i = \U(pn)\ 2 dp 1t . (4-4) 
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As fnipn), we take again a Gaussian function: 



Noting that 

d Pt ~ {E*/E°)dp« , Pi ~ p° + (p w - p») , 

where 

one may use Eq.(4.4) to translate f w (p w ) into fi(pi). This leads one to Eq.(2.8), with Api 
given by 

^ = (E°/E°„)Ap« , (4-6) 

which implies 

Api/Apj = E$/E$ , (4-7) 

to be compared with Eqs.(3.7). 

To appreciate difference between Eq.(3.7) and Eq.(4.7), it is instructive to note that, for 
the present case, the wave-packet factor ^(x,t), E.(2.13), may be expressed as 



W ^ t] -\^Ax^ [ 2(^^)2 



ex P[ 7^- A — W ] , O 8 ) 



where 



t!£\x,t) = m ± - (Pi ±P°)0K , (4-9) 



and /Aa^ = 1/Ap n , and that the velocity of the wave packets, or may be expressed 



as 



(Ef + E°y + (Ef - E»y ' 1 ' 



1 ,,0, _ (g + gXgg + p?) + (g - gXgg - p?) UU ) 

(rf+p°) 2 + (p°-p°) 2 " 



Since, for relativistic neutrinos, the second term is negligible compared to the first term in 
the numerator as well as in the denominator in each of Eqs.(4.10) and (4.11), we see that 
the center-of-mass velocity prescription is good enough to describe oscillations of relativistic 
neutrinos emitted from pions. 
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5. Discussion 



In Ref.l, we have, as another possible prescription, proposed the equal- velocity prescrip- 
tion and pointed out that this prescription would also lead one exactly to Eqs.(3.3) and 
(3.4)0and that, if the condition of equal- velocity indeed prevails, neutrino oscillations would 
be free from the coherence conditions, Eqs.(2.20) and (2.25). We have also mentioned that 
whether and how experiments to be described by the equal-energy/-momentum/-velocity 
prescriptions could become realistic and feasible remain to be carefully examined and con- 
trived with efforts. Here we like to mention that the equal-energy (-momentum) prescription 
seems appropriate to approximately describe oscillation experiments which involve something 
corresponding to time- (space-) integration or averaging.S 

For neutrinos produced in pion decays or in any two-body decays, none of equal-energy, 
equal-momentum and equal-velocity holds. ^Explicit kinematical considerations combined 
with our key postulate Eq.(4.4) lead us to Eq.(4.7) and hence to an explicit example which 
indicates that the center-of-mass velocity prescription is not necessarily applicable, or stated 
in a different way, the wave packets do not necessarily follow the classical trajectory described 
by Eq.(3.5). Although we have at the same time confirmed that this prescription is good 
enough to describe oscillations of relativistic neutrinos emitted from pions, possible deviation 
from a classical picture should be taken into account in general. 

To conclude, we like to remark that the main points discussed in the preceding sections 
and summarized above are approximation- dependent, and, for comparison, we shall give 
another way of approximation in Appendix B. 
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*) Note that /3° = (3% = p% = $ , if 0? = $ . 

**) The equal- velocity case was later considered also by De Leo et al..0 ) Ahluwalialilf suggested that the 
issue of equal- velocity versus equaLfinergy/-momentum can be quite important and should be testable for 
supernova neutrinos, while GiuntiElP suggested that, although the equal-energy/-momentum case could be 
realized approximately in some physical processes, the equal-velocity case seems very unlikely in any physical 
process. 

***) A similar remark was raised explicitly by Winter before. E3 
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Appendix A 

Remarks on normalization 

With Eq.(2.4), one has 

dx(u a/ (x,t)\u a (x,t)) = 5 a > a , (A-l) 



/°° / 9/W /\r • 

X 6XP| - 2 {( A.,f + ( A x ,f) 1 ' (A ' 2) 

/oo roo 
dxP a ^ oc '(xjt) = / dx{v a (X) t)\v a (X) t)) 
-oo J — oo 

= 1 . (A-3) 

Eq.(A.2) implies that, with our normalization (2.4), the initial condition stated in the be- 
ginning of Sect. 2 as "suppose a neutrino of flavor a is born at t = 0" is realized strictly 
only in the limit of Axj = Axj — > or if the momentum distribution functions fi(pi) have a 
common central value and a common dispersion. 

If, as a normalization condition alternative to Eq.(2.4), one adopts 




2vr 



dpifM = 1 , (A-4) 



one would have 



(v a ,(x,t)\v a (x,t)) = 5 a i a , (A- 5) 

F«-^(0,0) = (A-6) 

Y^Pa-*a'{x,t) = (V a (x,t)\v a (x,t)) 
a' 

= £tf«,^exp[-(s-/^)7(A*fl . (A7) 

i 

This way of normalization has advantage that it keeps the correspondence with the conven- 
tional plane- wave treatment as far and transparent as possible. In fact, this normalization 
ensures Eqs.(2.2) and (2.5) to reduce respectively to Eqs.(3.1) and (3.2), and, accordingly, 
J2 a ' Pa^a'(x,t) = (u a (x, t) \ v a (x, t) ) — > 1, when Api = l/Axi — > 0. Although the main points 
we discussed in the text are conceivably independent of how fi(pi) is normalized, the prob- 
lem of normalization itself involves some subtlety and ambiguity and deserves to be studied 
further.0 

*- > _Eq-(A.4) was adopted in our first paper, 0) and Eq.(2.4) has been adopted in our subsequent discus- 
sions. US' See Ref.7 for another normalization procedure. 
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Appendix B 

Another way of approximation 

To calculate Gij{x,t) (denned by Eq.(2.7)) for i ^ j, if one approximates the exponent 
(Pi - Pj)x - {Ei - E'j)t as 

(Pi - p'j)x - {Ei - E\)t 

= (ft - p'j)x - [(ft 2 + ™?) - (ft? + ™%)]t/(Ei + E' 3 ) 

= (Pi-P'j)[x - {p i +p' j )t/{E i + E' 3 )} - {ml - m 1 j )tl{E l + E[) 

~ { Pl - p}) [s - (p° + p°)*/(£° + - {m} - m))t/{E» + Efj) 

= {P. ~ P'j -P°i+ P?)(x - 0° t) + 0° (s, t) , (B-l) 

and performs integration, by changing the integration variables from pi and p'j to 

Pij = ft ~ P'j , Pa = {Pi^Y + p'^Ap^/^Apif + {A Pj f) , 

one would be led to 

a result compatible with the statement that the wave packets are on the classical trajectory, 
(3.5).0Note that Eq.(B.l) is to be compared with our approximation 

{ Pi - p'Jx - {Ei - E'j)t 

~ ( p . _ p '.) x _ [E o + ^ _ p oj _ E o _ pay _ p o )]t 

= {Pi j -P%)^-Plt)-{P ij -P^-^)t + 0%{x,t) , (B4) 
and Eq.(B.2) is to be compared with Eqs.(2.12) ~ (2.14), where 

P% = ft° ~ ft? , P% = (ft°(^ft) 2 + P^Ap^/^Apif + {A Pj f) , 
and ffij is defined by Eq.(2.16). 



*> In cotrast, in Ref.6, the authors presupposed that the neutrinos should be on their classical trajectory, 
and treated the phase facor (pi — p'j) — {Ei — E'A as 

{pi - p'^x - {Ei - fyt ~ -(m? - m 2 j )t/{Ef + E°) , (B-3) 

which is to be compared with Eq.(B.l). 
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